Nearly-logarithmic decay in the colloidal hard-sphere system 
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Nearly-logarithmic decay is identified in the data for the mean-squared displacement of the colloidal hard- 
sphere system at the liquid-glass transition [v. Megen et. al, Phys. Rev. E 58, 6073 (1998)]. The solutions of 
mode-coupling theory for the microscopic equations of motion fit the experimental data well. Based on these 
equations, the nearly-logarithmic decay is explained as the equivalent of a (3-peak phenomenon, a manifesta- 
tion of the critical relaxation when the coupling between of the probe variable and the density fluctuations is 
strong. In an asymptotic expansion, a Cole-Cole formula including corrections is derived from the microscopic 
equations of motion, which describes the experimental data for three decades in time. 

PACS numbers: 64.70.Pf, 61.20.Lc, 82.70.Dd 



The transition of a fluid to a glass is accompanied by a 
dramatic slowing down of the dynamics as the system ap- 
proaches the transition point |1]. The correlation function 
for some variable A, A (f) = (A* (t)A) / (\A\ 2 ) , with <> de- 
noting canonical averaging, decays rapidly typically within a 
time increase by a factor of 10 in the normal fluid regime, 
while it stretches over orders of magnitude in time close to 
a glass transition. Simultaneously nontrivial dynamical fea- 
tures like power law decay and unconventional scaling are ob- 
served. The characterization and explanation of these signa- 
tures of the glass transition present challenges to experiment, 
computer simulation and theory alike. One such feature is 
the nearly-logarithmic decay of the correlation function which 
was discovered recently in the orientational correlation func- 
tion by optical-Kerr-effect (OKE) measurements in molecular 
liquids [2]. The response function measured in these exper- 
iments is proportional to the negative time derivative of the 
correlation function, so a measured nearly-? -1 decay is equiv- 
alent to a nearly-logarithmic relaxation in 0a (0- m a subse- 
quent fit it was shown that the data can be welldescribed by a 
schematic model within the mode-coupling theory for ideal 
glass transitions (MCT) |3]. The nearly-logarithmic decay 
could be interpreted as a [3-peak phenomenon of the short- 
time critical relaxation when rotation-translation coupling is 
sufficiently large. Different from earlier suggestions, the pres- 
ence of higher-order glass transitions, like in systems with 
very short-ranged attractions 1 4], could be ruled out as a cause 
for the nearly-logarithmic decay in molecular systems. 

Schematic models are truncated versions of the full equa- 
tions of motion of MCT that share the mathematical structure 
and universal properties of the latter but do not contain micro- 
scopic details |5]. Schematic models involve fit parameters 
in the equations of motion for the correlation functions, what 
limits the interpretation of fits like in Ref. ||3|] with respect 
to the microscopic origin of the observed features - like the 
nearly-logarithmic decay - that go beyond universal predic- 
tions. In the full microscopic equations, these parameters are 
completely determined by the number density p = N/V for a 
system of /V particles in a volume V, and the static structure 
factor S q which is given by the interaction potential 10,0,0]. 
In the following, a system of hard spheres of diameter d shall 



be discussed, where the only external control parameter is 
the packing fraction cp = 7T,paf 3 /6. Predictions for the hard- 
sphere system (HSS) can be tested in colloidal suspensions, 
and the results for collective density fluctuations are reviewed 
in 1 8], demonstrating the universal laws of MCT and verify- 
ing certain characteristic parameters of the theory. The tagged 
particle dynamics has been analyzed recently in a computer- 
simulation study showing good agreement with MCT [9]. 

A particularly informative quantity to be studied in a system 
that shows structural relaxation is the mean-squared displace- 
ment (MSD), which is defined by 5r 2 (r) = (\r(t) - r(0)| 2 ) 
with r(t) denoting the position of a particle at time t. The 
MSD is studied frequently in computer simulation and was 
measured for nine decades in time in a colloidal suspension 
lUOll . In this paper, the MSD will be calculated for the HSS 
within MCT in order to show, that for a state near the glass 
transition there is a window of two decades in time of nearly- 
logarithmic relaxation. This will be done by identifying a f _1 
decay of the derivative of the MSD in complete analogy to the 
procedure used for the analysis of the corresponding result for 
the OKE data 1 2\. Second, by asymptotic solution of the MCT 
equations, a modified Cole-Cole law will be derived for the 
description of the critical decay of the MSD. This analytical 
formula explains the nearly-logarithmic regime and accounts 
for the critical decay for three decades in time. Third, it will 
be shown that the MCT solutions at the critical point describe 
the experimental result in 1 10] for five orders of magnitude in 
time. Fourth, the full range of available experimental data for 
the MSD is described reasonably by the numerical solutions 
of the MCT equations going beyond the universal laws. Thus 
nearly-logarithmic decay laws known from molecular systems 
as p-peak phenomenon are established in a colloidal system. 

Slow dynamics within MCT originates from a singularity 
in the long-time behavior of the collective density fluctuations 
|0]. Other variables share that behavior if they are coupled to 
the density fluctuations. While the coupling to rotational mo- 
tion was modeled by an empirical parameter in |3], the MSD 
for a colloidal system is given be a microscopic equation of 
motion by 
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FIG. 1 : Time derivative of the mean-squared displacement D(t) in 
the hard-sphere system for packing fractions cp = cp e = 0.5159 and 
cp = 0.5145 (full lines from left to right). The upper dashed line is 
proportional to t . The lower dashed curve represents the approx- 
imation by Eq. 0. The dotted line displays the critical law t~" . 



Here and in Fig.|2|the unit of time is fixed by the short-time diffusion 
coefficient D^/d 2 = 1/160. The inset shows D(t) for cp = cp c on a 
linear scale for comparison. 

with the short- time diffusion coefficient D s {) . The memory ker- 
nel m^\t) is determined by the static structure of the liq- 
uid and collective and single density correlation functions. 
The structure factor is calculated using the Percus-Yevick 
approximation (PYA) |6], details of the numerical solution 
are outlined elsewhere II ill . For cp below the critical point 
cp c = 0.5159, the long-time solution of Eq. ( fTal is diffusive, 
5r 2 (t) = 6D s t, with a long-time diffusion coefficient D s de- 
pending on the distance from the critical point. For cp ^ (p e , 
the long-time solution is arrested at a plateau 8r 2 (f) = 6r 2 , 
with a localization length r s . The derivative of the MSD 
can be interpreted as a time-dependent diffusion coefficient, 
6D(t) = d[5r 2 (f) / d 2 } / dt , that goes to zero as the glass tran- 
sition is approached \12]. D(t) calculated from Eq. dTab is 
plotted for <p = 0.5145 and <p = <p c as full lines in Fig.^ The 
dynamics when approaching the critical point shares more 
and more of the relaxation at cp c as seen for cp = 0.5145 for 
log 10 f < 1. This portion of the dynamics is therefore referred 
to as critical relaxation. For times exceeding the critical re- 
laxation, the dynamics crosses over to the long-time diffusion 
as seen for (p = 0.5145 around log 10 f w 2. The most signif- 
icant finding in connection with Fig. \l\ is the appearance of 
a window in time after the transient where the dynamics fol- 
lows closely a t law for logi f ^ 2. A comparison of 
the solution for cp = 0.5145 with the critical decay shows that 
the major fraction of the f _1 decay is part of the critical re- 
laxation. This decay reflects the same behavior as found in 
molecular liquids, moreover, the decay observed in Fig.^is 
closer to t than for most of the OKE data where t ~ x is seen 
with x ranging from 0.8 to 1.15 The long-time decay at 

the critical point is a r _a -law with a=0.312 for the HSS. This 
law is shown as dotted straight line of slope — (1 +a). It ac- 
counts for the critical decay only for log 10 f 3, is preceeded 
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FIG. 2: Mean-squared displacement at the critical point in the hard- 
sphere system for Newtonian (N) and Brownian (B) dynamics (full 
curves). The short horizontal line displays the plateau 6r£ 2 . Dashed 
lines show the asymptotic approximations by Eqs. J2j (labeled t~ 2a ) 
and J4}. The leading-order approximations are shown dotted for 
Eq. ^2) (t~ a ) and Eq.^4j (cc). The nearly-logarithmic relaxation is 
shown by the straight dashed line labeled In?. The diamond and the 
square indicate 10% deviations of the critical solution B from the ap- 
proximation by Eq. and the solution for cp = 0.5 145, respectively. 

by the f~' decay, and has no relevance for <p = 0.5145. 

The nearly-logarithmic decay is found adjacent to the tran- 
sient for the Brownian colloidal system as for the Newtonian 
molecular systems [3]. For the HSS also Newtonian dynamics 
can be considered and the equivalent of Eq. ( fTal reads 1 13], 

8 r 5r 2 (f) + v 2 f dt'm^{t-t')br 2 {t')=6v 2 s t, (lb) 
Jo 

with the thermal velocity v 2 for a tagged particle. Figure |2 
shows the solution for the MSD for Newtonian and Brownian 
dynamics for cp = cp 6 . Both solutions can be matched at long 
times and they overlap down to 8r 2 (f) w 0.01 defining the be- 
ginning of the regime for structural relaxation at log 10 f w 
11311 . The nearly-logarithmic decay is shown as nearly-linear 
increase of the MSD with logf outside the transient within the 
same time interval as marked by the dashed line in Fig^ 

To further analyze the window of nearly-logarithmic relax- 
ation and in an attempt to distinguish it from a further approx- 
imation to a power law with very small exponent like in 1 14], 
asymptotic expansions shall be applied. At the critical point, 
the MSD can be expanded in power laws and reads ITUl . 

8r 2 (0/6 = r c2 - h MSD (t/t )- a + liMSDKMSD{t/tor 2a , (2) 

where the first two terms on the RHS constitute the leading- 
order result. The values for the HSS are r c H = 0.0746, Iimsd — 
0.0116, K MS D = -1-23, t = 0.425, X = 0.735. Leading 
and next-to-leading order result of Eq. (|2ji are shown dotted 
and dashed in Fig. [2] and describe the full solutions only for 
t > 1000 and t > 100, respectively. As will be shown below, 
the leading order fails to apply in the experimentally relevant 
window, while including the correction still does not explain 
more than one decade of structural relaxation. 
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In the following, a new asymptotic solution is derived 
which is based on the expansion of the memory kernel. For 
long times, the equations of motion for the MSD are the 
same for both Newtonian and Brownian dynamics ll3Tl . 
and can be represented with the modified Laplace transform 
S \f(t)](z) = izJodt exp[/zf]/(f) as 



5[5r 2 (f)](z) = 



5[m(°)(f)](z) 



(3) 



An asymptotic expansion for the memory kernel equivalent to 

the one in Eq. (0 is given by m( 0) (f) = / c (0) +h m{0) (?/?o)~ a + 

('Ao) -2 ". with f m {o) = lAf, h m(0) =h MS D/rf, 
K m (o) = Kmsd + UimsdK 2 , and X = T(l - a) 2 /T(l - 2a) 
with the Euler Gamma function T(x). In the Laplace do- 
main, this can be written with a characteristic frequency Cflp as 

S[m(°\t)}(z) = /; n(0) {l + (-/z/cop) u + ^« SD (-/z/co p ) 2fl }, 
where K™ SD = f c m(0) K m(0j /(h ml0) X). Hence, Eq. Q reads up 
to next-to-leading order 



S[5r 2 (t)](z) = 



6r! 



c2 



i + (-iz/^y+K^ SD (-iz/®v) 2 « 



(4) 



and 

to 



frequency is given by cop 



the characteristic 

j^p^ lA '. For the HSS we get r cop = 0.03895. The 
leading-order result in Eq. ® is obtained for K^ SD = 0, and 
is known as Cole-Cole law. For different variables a similar 
result was obtained before SHI. While the correction in 
Eq. (0 is comparably large, Kmsd = — 1 -23, the analogous re- 
sult for the memory kernel is only a fraction of it, K ' = 0.30. 
The correction in Eq. (0} is even smaller, K C ^ SD = 0.196. As 
the corrections determine the range of validity for the leading 
terms of the asymptotic expansions, the leading-order result of 
Eq. @ is superior to the leading-order result of Eq. (0 as seen 
in Fig. where already the leading order of Eq. (0} describes 
qualitatively the complete relaxation. Including the correction 
explains the complete window of structural relaxation as seen 
by the diamond marking a 10% deviation of the approxima- 
tion by Eq. (0J from the solution of Eq. ( fTal . A comparison in 
Fig- n shows that Eq. (|4} also covers the complete regime of 
nearly-logarithmic decay. For (p = 0.5145, D(t) in Fig.[TJand 
the MSD in Figgis described by Eq. gj for log 10 f < 3. The 
crossover from the critical decay to the long-time diffusion 
yields an extended window in time where lnf fits the solution. 

Figure [3] shows the data from 1 1011 in units of the particle 
radius R = d/2 and the Brownian time scale T& = R 2 / (6Z)q) 
together with the solutions of Eq. ( fTal for different values of 
(p. The agreement is rather satisfactory considering that all 
coupling parameters are fixed microscopically. Three adjust- 
ments need to be made to match experimental and theoreti- 
cal curves. First, the short-time diffusion coefficient sets an 
overall time scale that needs to be accounted for. The the- 
oretical solutions are matched to the data for log 10 f < — 1. 
As the experimental curves do not fall on top of each other 
in that regime, the chosen time scale also needs to vary by a 
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FIG. 3: Fit of the mean-squared-displacement data from [J^J (di- 
amonds and full curves) by the solutions of mode-coupling theory 
(dashed). The dotted curve displays the theoretical result at the crit- 
ical point and is equivalent to curve B in Fig. [2] The square indi- 
cates a 10% deviation between the critical relaxation and the fit for 
(p ex P = 0.583. The open diamond is the same as in Fig.|2| The in- 
set shows the mapping of the experimental packing fractions (p exp 
onto theoretical values (p theo on the left axis as closed and open cir- 
cles for curves shown or left out, respectively. The dotted line shows 
the function (f ±eo = 0.87 (p ex P + 0.023. The right axis displays the 
vertical scaling factor by crosses. 

factor up to 2. Second, with cp 6 ' « 0.58 CUE], the experi- 
mental glass transition in the HSS is found at approximately 
12% higher packing fraction than predicted in Refs. 171 ll ill . 
This discrepancy could be reduced to below 10% by using the 
actual structure factor rather than the PYA in a computer sim- 
ulation 1 17], but a remaining error from the mode-coupling 
approximation has to be accepted. This is similar to other ap- 
proaches fitting the data 1 1 8 ] . Therefore, theoretical values for 
(p are mapped on experimental ones by comparing the diffu- 
sivity in experiment and theory. The result of this procedure 
is shown in the inset of Fig.[3]by circles. The mapping is fitted 
by (p theo = 0.87<p exp + 0.023, which is remarkably similar to 
what was found from a computer simulation of the HSS j^l, 
and extrapolates to a critical density <p 6 « 0.57. The highest 
values of cp do not follow that linear mapping due to ageing 
for log 10 f/T£ > 4 fiolflrjl . The dynamics before ageing sets 
in, and hence the regime for nearly-logarithmic relaxation, is 
not affected significantly by ageing 1 16]. Third, as the actual 
transition takes place at a higher packing fraction than the- 
ory predicts, a smaller localization length and hence a vertical 
shift of the curves by a prefactor a y - shown by crosses in the 
inset of Fig. [3]- can be expected. Surprisingly, many curves 
are described well using a y = 1 . The variation of a y yields an 
error estimate for the localization length of around 20% for 
the localization length as compared to the theoretical results. 
If only data close to the glass transition is considered, the data 
show indeed a smaller localization length than predicted. Af- 
ter these adjustments, Eq. ( flat fits the data for different pack- 
ing fractions over the entire experimental range of up to nine 
orders of magnitude in time. 



4 



The curve (p ex P = 0.583 is fitted by (p theo = 0.5146, and in 
addition, the critical relaxation for (p theo = cp e is matched with 
the same prefactors and shown dotted. The data follow closely 
the critical relaxation for t/z/, < 2000. The two theoretical 
curves correspond to the ones shown in Figs. ^ and [2] thus 
establishing the existence of a window of nearly-logarithmic 
relaxation in the experimental data. 

In conclusion, it is shown in Fig.[5]that mode-coupling the- 
ory (MCT) is able to describe the mean-squared displacement 
(MSD) in a colloidal hard-sphere system (HSS) 1 10] rather ac- 
curately. This implies consistency with both universal laws of 
MCT and non-universal MCT predictions for the HSS. Out- 
side the transient, the data exhibit a window of two orders 
of magnitude in time where the relaxation resembles nearly- 
logarithmic decay or a p-peak phenomenon, cf. Fig. ^ The 
latter is found independent of the underlying Brownian dy- 
namics also for Newtonian dynamics, cf. Fig. [2] and is distin- 
guished from an accidental crossover by an asymptotic expan- 
sion relating it to a p-peak, cf. Eq. ©. Such a dynamical fea- 
ture was discovered in molecular systems |2] and interpreted 
as a consequence of strong translation-rotation coupling |3]. 
In addition, similar nearly-logarithmic decays were also ob- 
served in nematic liquid crystals |19jy. For the hard-sphere 
system (HSS), while clearly lacking rotational degrees of free- 
dom, the MSD constitutes another variable that is coupled 
strongly to the density fluctuations if the tagged particle is suf- 
ficiently large. It can be inferred from Ref. II ill that m(°> (t) in 
Eq. Q becomes smaller when the size of the tagged particle 
is chosen smaller than the particles in the host fluid. A larger 
particle in turn experiences a stronger coupling of its MSD to 
the collective dynamics. Hence, the counterintuitive finding in 
01, that the dynamics of a coupled variable can deviate further 
from the underlying collective dynamics as the coupling is in- 
creased, also applies to the MSD of the HSS for increasing 
the diameter of the particle for which the MSD is measured. 
The empirical nearly-logarithmic law is part of the critical re- 
laxation and does not change slope when observed sufficiently 
close to the transition point. This is in clear contrast to the log- 
arithmic laws in the vicinity of higher-order glass-transition 
singularities 14] and can be used to distinguish both relaxation 
features. Hydrodynamic interactions present in colloidal sus- 
pensions are not included in Eq. {B but are known to change 
the short-time dynamics in MCT 12011 . The quality of the data 
fit and the validity of the nearly-logarithmic law are therefore 
surprising and rule out a significant contribution of hydrody- 
namic interactions on the dynamics of the MSD within MCT 
as shown in Fig. which is in contrast to other approaches 
This can be rationalized by the strong coupling of 
the MSD to the collective dynamics causing the MSD for a 
larger particle to be slower than for a smaller particle 11 ill . 
For the slowest collective relaxations only small deviations by 
hydrodynamic interactions have been found 1 20], a result that 
apparently seems to apply also to the slow coupled variable 
MSD. 

Equation @ improves the understanding of glassy dynam- 
ics in three aspects. First, it comprises a Cole-Cole that is de- 
rived rigorously from a microscopic equations of motion and 



therefore provides a foundation for data interpretations as in 
Ref. Q]. Second, it closes a gap in the understanding of the 
MSD close to the glass transition 1 1 311 . Together with the re- 
sults from Ref. 1 1 1], MCT provides analytical descriptions of 
the dynamics in all subsequent windows in time up to small 
crossover regions: The initial ballistic or Brownian dynamics 
is followed by a critical relaxation described by Eq. @, before 
a plateau 6r^ 2 is reached. After that plateau, the crossover 
to the long-time diffusion starts with a von Schweidler law 
r b where b = 0.583 Q. Third, precise meaning is given to 
the finding that strongly coupled variables do not necessarily 
show the same dynamics. 
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